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We calculate the transverse momentum Q± dependence of the helicity structure functions for the 
hadroproduction of a massive pair of leptons with pair invariant mass Q. These structure functions 
determine the angular distribution of the leptons in the pair rest frame. Unphysical behavior in 
the region Q± —* is seen in the results of calculations done at fixed-order in QCD perturbation 
theory. We use current conservation to demonstrate that the unphysical inverse-power and ln(Q /Q±) 
logarithmic divergences in three of the four independent helicity structure functions share the same 
origin as the divergent terms in fixed-order calculations of the angular-integrated cross section. We 
show that the resummation of these divergences to all orders in the strong coupling strength a s can 
be reduced to the solved problem of the resummation of the divergences in the angular-integrated 
cross section, resulting in well-behaved predictions in the small Q± region. Among other results, 
we show the resummed part of the helicity structure functions preserves the Lam- Tung relation 
between the longitudinal and double spin-flip structure functions as a function of Q± to all orders 
in a s . 



I. INTRODUCTION 

Production of a massive pair of leptons of opposite elec- 
tric charge in hadronic interactions, hi + hi — > l + (r X , 
has revealed new narrow hadronic states, notably the 
J and the T, and it continues to provide an important 
complement to deep-inelastic lepton scattering and other 
hard-scattering processes for probing the short-distance 
dynamics of strong and electroweak interactions. The as- 
sumption that the broad continuum of £ + £~ pairs origi- 
nates from quark- antiquark annihilation through a single 
virtual photon, as embodied in the Drell-Yan model [lj, 
implies that the angular distribution in the £ + £~ rest 
frame should be that of a transversely polarized photon, 
(1 + cos 2 8), where the polar angle 9 is the direction of the 
lepton relative to the direction of the incident quark and 
antiquark. Acceptance restrictions limit measurements 
of the full angular distribution, but qualitative verifica- 
tion of this expectation was one of the early tests that 
increased confidence in the model [2|. 

In practice, massive lepton pairs are produced with 
substantial transverse momentum supplied from 

a theoretical perspective by higher-order processes in 
perturbative quantum chromodynamics (QCD). An in- 
teresting challenge has been to predict how the an- 
gular distribution should behave as a function of Q± 
E II H & S B 0, EE EH- Indeed, this challenge is 
part of the more general ambition to predict the fully 
differential cross section da/dQdQ±dyd$l, where Q is the 
invariant mass of the lepton pair, y is its rapidity, and 
dfl = d cos 8d<fi represents the differential decay angular 



distribution in the pair rest frame with respect to a suit- 
ably chosen set of axes. 

In addition to the virtual photon, the W boson and 
the Z boson also have important decay modes into pairs 
of leptons. The angular distribution of these leptons, 
measured in the rest-frame of the parent states, deter- 
mines the alignment (polarization) of the vector boson 
and, consequently, supplies more precise information on 
the production dynamics than is accessible from the spin- 
averaged rate alone. An understanding of the changes 
expected in the angular distribution as a function of the 
transverse momentum Q± is a topic of considerable im- 
portance, both for refined tests of QCD and for elec- 
troweak precision measurements. An example of a QCD 
process is the flavor dependence of W production in po- 
larized hadron-hadron scattering at the Brookhaven Rel- 
ativistic Heavy Ion Collider (RHIC) [Hj]. Better under- 
standing of the expected angular distributions will reduce 
the systematic uncertainties on the determination of the 
W boson mass [l3l [lij and, in turn, improve the bound 
on the mass of the Higgs boson within the standard model 
of particle physics. 

In this paper we consider the scattering of two hadrons 
of momentum P\ and P2, respectively, producing a vir- 
tual photon of four-momentum q, A{P\) + B(Pi) — > 
r y*(q)+X, that decays into a pair of leptons of momentum 
I and I, as sketched in Fig. [T] The ideas and techniques 
developed here can be applied readily to the production 
of W and Z bosons, as well as to other yet-to-be-observed 
massive vector bosons that decay into a pair of leptons. 
They are applicable also in semi-inclusive deep-inelastic 
scattering (SIDIS). 
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The general formalism for the description of the an- 
gular distribution in terms of helicity structure functions 
is developed for the Drell-Yan process in Ref. [4|. The 
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differential cross section may be expressed as 



da 



d 4 q dn 



[w T {i 



2(2tt) 4 S 2 Q 2 
+W L (l-cos 2 6) 



W A (sin(26>) cos. 



-W A A(sin 2 6>cos(20))] 



(1) 



The four independent "helicity" structure functions Wt, 
Wl, Wa, and Waa depend on Q, Q±, rapidity y, and 
on the center-of-mass energy \fS of the production pro- 
cess. They are defined in the virtual photon's rest frame, 
and they correspond, respectively, to the transverse spin, 
longitudinal spin, single spin-flip, and double spin-flip 
contributions to the Drell-Yan cross section. 

The angular-integrated cross section is expressed in 
terms of Wt and Wl as 



or 



da 



d 4 q 12ir 3 S 2 Q 2 



[2W T + W L ] 



(2) 



An interesting relationship Wl = 2 Waa between the 
longitudinal and double-flip structure functions is derived 
in Ref. [4j in the context of the parton model, and it 
has been shown to hold at least approximately at higher 
orders in perturbative QCD. Ex per imental tests of this 
relationship are reported in Ref. jl5l . H5 , Il7l . [l8l [T^ . 

Our principal focus in this paper is the prediction of 
the full Q± dependence of the four structure functions, 
including the region of small and intermediate Q± where 
the cross section takes on its largest values. Many papers 
dealing with various aspects of Drell-Yan angular distri- 
butions have preceded ours. Explicit perturbative calcu- 
lations were done in the parton model 0,0, in perturba- 
tive QCD at order a s pi, [g, M , 3 and a 2 0, as well as in 
high twist formalisms [20 . [21| . When calculated at fixed 
order in QCD perturbation theory, the structure func- 
tions show unphysical inverse-power QJ™ (n = 1 or 2) or 
logarithmic ln(Q/Q±) divergences, or both, as Q± — ► 0. 
For the angular-integrated cross section, da/d 4 q, it is well 
established that similar unphysical divergences can be re- 
moved after resummation of the ln m (Q 2 /(5 2 L )/(5 2 L singu- 
lar terms from initial-state gluon emission to all orders 
in a s HI HHHHI. 

Examinations of the singular logarithmic terms in the 
helicity structure functions are reported in Refs. [l(J HH, 
[11, Eg. Since only W T shows the ln m (Q 2 /Q 2 L )/Q 2 L di- 
vergence, previous resummation calculations were carried 
out only for Wt in the same way as for the angular- 
integrated cross section. As shown in Refs. [H, Ha, 
resummation removes the perturbative power divergence 
in Wt- One consequence of resummation of just Wt is 
a large change in the relative size of Wt and the helic- 
ity structure functions for which no resummation is per- 
formed. This result is not quite consistent with general 
expectations about the relative size of helicity structure 
functions in the Collins-Soper frame. For example, one 
expects Waa/Wt -> Q\ as Q± -> HJ. 

In Ref. [ll|, Boer and Vogelsang carefully investigate 
the logarithmic behavior of the order a s perturbative 



contributions to the helicity structure functions. At or- 
der a s , they find that, like Wt, both Wl and Waa 
have a ln(Q 2 /Q 2 L ) logarithmic divergence, but not the 
power divergence seen in Wt, and that Wa has no 
logarithmic divergence at this order in the Collins-Soper 
frame. They notice that the logarithmic contribution to 
Wl and Waa from quark-gluon (or gluon-quark) sub- 
process is different from that for Wt and does not fit the 
pattern expected for the perturbative expansion of the 
Collins-Soper-Sterman resummation formalism to order 
a s [25|. They also discuss the frame dependence of this 
logarithmic contribution. 

The present paper expands on our earlier short 
manuscript on the same subject We start with 

the observations that the four helicity structure func- 
tions cannot be independent at Q± = and that the 
general tensor decomposition in the virtual photon rest 
frame in Eq. ^ is ill-defined at Q± = 0. Then, guided 
by electromagnetic current conservation, we construct a 
new asymptotic form for the hadronic tensor with the 
right degrees of freedom as Q± — > 0. We find that 
the leading logarithmic behavior of the different helicity 
structure functions, Wt, Wl, and Waa, has a unique 
origin. We reduce the problem of transverse momen- 
tum resummation for Wt, Wl, and Waa to the known 
solution of transverse momentum resummation for the 



angular-integrated cross section |25| , and we prove that 
the logarithmic divergences in Wt, Wl, and Waa may 
be resummed to all orders in the strong coupling strength 
a s , yielding well behaved predictions that satisfy the ex- 
pected kinematic constraints at small Q± . We emphasize 
three main results of our research: 

• Current conservation uniquely ties the perturbative 
divergences as Q±/Q — > of the otherwise inde- 
pendent helicity structure functions Wt,Wl, and 
Waa to the divergence of the angular-integrated 
cross section. 

• The perturbative divergence in the angular- 
integrated cross section is sufficient to remove all 
leading divergences of the four individual helicity 
structure functions. 

• Transverse momentum resummation of the 
angular-integrated cross section determines the 
resummation of the large logarithmic terms of the 
helicity structure functions Wt, Wl, and Waa, 
and the approximate Lam- Tung relation is an 
all-orders consequence of current conservation for 
the leading perturbatively divergent terms. 

The rest of our paper is organized as follows. In Sec. II, 
we define the helicity structure functions, and we derive 
the QCD perturbative contributions at order a s . We 
work in this paper entirely in the context of collinear 
QCD factorization [281], meaning that Q± > Aqcd, al- 
though Q±/Q may be small. We examine in detail 
the leading behavior of the perturbative contributions 
to the helicity structure functions in the limit of small 
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FIG. 1: Diagrammatic representation of hadronic dilepton 
production via a virtual photon of four-momentum q. 



Q±/Q. In Sees. Ill and IV, we investigate the generic 
singular structure of the perturbative contribution to the 
Drell-Yan hadronic tensor, and we derive an asymptotic 
current-conserving tensor that explicitly includes all the 
leading divergences of the perturbatively calculated he- 
licity structure functions in the limit Q±/Q — > 0. We 
also explore the connection between cross sections for in- 
cident parton states of fixed helicity and the subleading 
perturbative contribution to the spin-averaged helicity 
structure functions. We discuss all-orders transverse mo- 
mentum resummation for helicity structure functions in 
Sec. V obtaining well-behaved distributions as a func- 
tion of Q±. We show that the resummed part of the 
helicity structure functions satisfies the Lam- Tung rela- 
tion, Wl = 2WaAj between the longitudinal and the 
double-spin- flip structure function to all orders in a s . 
Finally, in Sec. VI, we offer a summary and our con- 
clusions, and we outline plans for future work on W and 
Z hadroproduction and in semi-inclusive deep-inelastic 
scattering. Three appendices are included in which we 
present detailed technical derivations of points discussed 
in the main body of the text. 



II. HELICITY STRUCTURE FUNCTIONS AND 
PERTURBATIVE CONTRIBUTIONS 

In this section, we define the helicity structure func- 
tions of Eq. ([T]). We present the next-to- leading order 
perturbative contributions to these functions and exam- 
ine the structure of the singular behavior of each helicity 
structure function Wj as Q±/Q — ► 0. 



Definition and normalization 



Helicity structure functions are defined in the virtual 
photon's rest frame. Let e^(<?) be the virtual photon's 
polarization vector with three polarization states, A = 
±1,0. The helicity structure functions are 



W T 






W l 


= Wp, 




W a 


= IV^(« + «)/V2, 




Waa 




(3) 



for the transverse spin, longitudinal spin, single spin-flip, 
and double spin-flip contributions to the Drell-Yan cross 
section, respectively. In the virtual photon rest frame 
(the center-of-mass frame of the dilepton pair), the po- 
larization vectors can be expressed in terms of orthog- 
onal unit vectors in that frame, V M , and Z^, as 
4 = (TV^ 1 - iyy-v/2, 4 = Z v fl. These unit vectors 
are normalized as X 2 = Y 2 = Z 2 = — 1, and they are 
also orthogonal to the current vector q^. They conserve 
the current, q^X^ = q^Y^ = q^Z^ = 0. Naturally, we 
can choose the fourth unit vector for the q = Lorentz 
frame to be T M = qf/Q with T 2 = 1 and Q = ^/q 2 . The 
full Drell-Yan hadronic tensor can be written in terms of 
the helicity structure functions and unit vectors in the 
virtual photon rest frame as [H 



T»T V ) {W T + Waa) 
-2X"X v Waa + Z"Z U {W L - W T 
- {X»Z V + X v Z^)Wa ■ 



Waa) 



(4) 



Different choices of the axes lead to different q = 
frames We choose to work in the Collins-Soper frame 
12611 . whose unit vectors are defined as, 



Z>' 



q P2 Pf - q Pl P% 



VQ 2 + Ql 

Y v = e ^ aP T„Z a X f3l 



qp 2 Pf + q Pl P* 



(5) 



where the dimensionlcss current conserving hadron mo- 
menta are P> = [Pf - (P • q)/q 2 q^/y/S with i = 1, 2, 
and qp { = Pi-qj \[S with i = 1,2. We present our deriva- 
tion and predictions on helicity structure functions in this 
Collins-Soper frame. Transformation of our results to 
other commonly used frames is simply a rotation around 
the V-axis H[lH. 

When the virtual photon mass Q and its transverse 
momentum Q± are much larger than Aqcd, we expect 
QCD collinear factorization to be valid for the Drell-Yan 
cross section (28[. Neglecting the transverse momentum 
of partons participating in the hard collisions, we write 
the incident parton momenta as 



a 
Pi 



LI 

P<2 



6P 2 M - 



(6) 



Neglecting all corrections suppressed by powers of 
Vqcd/Q or Aqcd/Q_l, we can factor the hadronic tensor 
as 



E 



6 



x< 6 ^ 7 »x(6,6,g), 



(7) 



with incoming parton distributions </>/(£) of flavor / and 
momentum fraction £. The short-distance partonic ten- 
sor is 



1' 



b->y*X 



x(2ir)U 4 { Pl +p 2 -q- 
d 3 p x 



n 



(8) 



At the most basic level, a massive virtual photon arises 
from quark-antiquark annihilation q + q — > 7* in a col- 
lision of hadrons, and it is produced with Q± = 0. The 
corresponding partonic tensor is 



1 



iCn" -I- n^ri" 



X£l6£(£l -^2) 

x (2ir) 4 S5 2 (Q ± ), (9) 

with color factor 1/3 and fractional quark charge e q . The 
unit vectors are n M = <5 M+ and rt^ 1 = and 



£1 = ~i= c 



X2 = 



(10) 



The lowest order helicity structure functions from qq 

7* are 



W, 



(0) 



< ) = <i = o. 



(11) 



First-order gluon radiation supplies finite through 
the quark-antiquark and quark-gluon subprocesses, q + 
q — > 7* + <? and g + g — > 7* + q, as sketched in Figs. [Hand 
|3l respectively. Perturbatively, these finite-order subpro- 
cesses yield singular differential cross sections as a func- 
tion of Q± in the limit Q±/Q — > 0. For the angular- 
integrated cross section, d<j/d 4 q, it is well established 
that this unphysical divergence can be removed after re- 
summation of the singular terms from initial-state gluon 
emission to all orders in a s [2||. The dependence of the 
helicity structure functions on Q j_ is our central focus in 
the rest of this manuscript. 



B. Order a 3 contribution 

In this section we present explicit expressions for the 
contributions at order a s to the four helicity structure 
functions from the two subprocesses qq — > 7*17 and 
qg — ► 7*q. Although some of the perturbative results 
are available in the literature, we present for complete- 
ness in Appendix IB1 and [C| the details of the perturbative 
calculation in a consistent notation for the spin-averaged 
and "polarized" contributions to the parton-level helicity 
structure functions in the Collins-Soper frame. 

To better identify the analytic behavior as Q±/Q — > 
of the perturbative contributions, we express the results 
in terms of two new variables, 



JZl = 



Z 2 



X2 

6 



(12) 
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FIG. 2: Feynman diagrams for quark-antiquark annihilation 
to a virtual photon plus a gluon. 



The parton-level Mandelstam variables defined in 
Eq. (|B3|) in Appendix [Bl are expressed as 



zi z 2 



Q 2 



1- Wi + QW 

i-z lv /i + Qi/Q 2 



(13) 



and 



1 

/; u 
1 



1 



S(-t) 



1 

iQl 



1- Wl + Qi/Q 2 



In the following subsections, we present our calculation 
for spin- averaged and polarized incident partons, with 
our specification of polarized states presented below. 



1. Spin averaged quark-antiquark annihilation 

As derived in Eq. (|B5[) in Appendix[Bj the contribution 
to the parton-level helicity structure functions from the 
quark-antiquark annihilation subprocess, after averaging 
over the spins of the incident quark and antiquark, are 



3 



s 

x S(s + t + u - Q 2 ) . 

ziz 2 



W L = e q —3— C F [Zi + Z 2 



s 

x 5(s + i + u-Q 2 ). 

Z\Z 2 



qq 2 87T 2 q s 1 

'■>AA " e q 3 2 



C F [z 2 + z 2 ] 



S 

x 5(s +i+u- Q 2 ) 

Z\Z 2 



■ w 



qq 

L ' 



_ >2 Sr-n. /' Q 



3 VQ 
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FIG. 3: Feynman diagrams for quark-gluon scattering to pro- 
duce a virtual photon plus a quark. 



S 

Z\Z 2 



5(s + t + u-Q 2 ). 



(14) 



The color factor is written as 4/9 = (1/3) x Cf, with 
(1/3) being the color factor for the lowest order con- 
tribution in Eq. ©, and C F = 4/3. With the ex- 
change of z\ and z 2 (or t and u), Eq. (fT4")) is also valid 
for the antiquark-quark scattering subprocess, except for 
w A q which acquires an extra overall minus sign that 
arises from the minus sign in the expression for wa in 

Eq. Kxm . 

The phase space (5-function can also be expressed in 
terms of the new variables as 



S 

Z\Z 2 



5(s + t + u-Q 2 ) 



1 



X\X 2 



S[ (l-z lx /l + Qi/Q 2 ) 



x(l - ^l + Qi/Q 2 ) - Qi/s) . (15) 



2. Spin-averaged quark gluon scattering 

As derived in Eq. (|B8[) in Appendix [Bl the contribu- 
tions from the quark-gluon subprocess, after an average 
over the spins of the initial quark and gluon, are 



,,qg _ 



2 87r 2 a s fQ^_ 

Ql 



l- Wi 



i 2 /Q 2 



[z 2 + ( Zl z 2 - If 



S 



Z\Z 2 



<5(s + i + -u-Q 2 ), 



,,95 _ 



87r 2 a s 



3 

x T R [z 2 + ( Zl + z 2 ) 2 ] 



S 



Z\Z 2 



5(s + i+u-Q 2 ). 



,,19 _ 



87r 2 a s 



i-Wi + Qi/Q 2 



x-T R [4 + {z x + z 2 ) 2 ] 

s 

x 5(s + i + u-Q 2 ) 

z\z 2 



1 



w ■ 



e 



2 87r 2 a s / Q 



1-Wl + QVQ 5 



x Tr [z 2 - 2z 2 2 ] 
x S(s + i+u-Q 2 ). 



Z\Z 2 



(16) 



The color factor is written as 1/6 = (1/3) x Tr with 
T R = 1/2. With zi and z 2 switched, Eq. (fl~6|) is also 
true for the gluon-quark scattering subprocesses, except 
for w A q which acquires an extra overall minus sign that 
arises from the minus sign in the expression for wa in 
Eq. (IAT21) . 



3. Expressions for polarized incident partons 

The behavior at small Q±/Q of the parton-level helic- 
ity structure functions is sensitive to the helicity states of 
the incoming partons. We present here the perturbative 
contribution to the helicity structure functions from the 
q+q — ► 7* + <? and q+g — ► 7* +q subprocesses with initial- 
state (anti)quark and gluon in a fixed helicity state. For 
a quark of momentum p, the helicity projection operator 



P±(p) = 7}1-P± t;1-PJ5, 



(17) 



where the first term on the right-hand-side (RHS) corre- 
sponds to the projection for a spin-averaged quark state, 
while the second term corresponds to the projection for a 
polarized quark state, defined as the state with incoming 
quark polarization projected onto the difference of the 
quark's helicity states. Similarly, the helicity projection 
operator for a massless gluon of momentum p moving in 
either the light-cone "+" or "— " direction is 



1 



1 



Pf(p) = -d a ? ±-ie af ' 



(18) 



where the transverse tensor d a ^ = —g a/3 + n a rfi + n a n^ , 
£ q/3 _ e a Pp a f lpUa ^ th e nrs t term on the RHS again cor- 
responds to the projection for a spin-averaged and phys- 
ically polarized gluon state, while the second term cor- 
responds to the projection to a polarized gluon state, 
defined as the state with incoming gluon polarization 
projected onto the difference of the gluon's physically 
polarized states. 

The contribution with a mixed unpolarized and a po- 
larized parton state leads to an antisymmetric contribu- 
tion to the hadronic tensor W^vi and it does not con- 
tribute to the Drell-Yan angular distribution. The sum or 
difference of our unpolarized and polarized contribution 
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correspond to the contributions from initial-state partons 
of the same or different fixed helicity state. 

Equation (|C1|) in Appendix [Cl shows that the polarized 
quark-antiquark contributions to the helicity structure 
functions are the same as the unpolarized contributions, 



u T , 



J L 1 
AA 



(19) 



with all unpolarized contributions given in Eq. |T] 

In treating quark-gluon scattering, we present results 
separately for the quark gluon and gluon quark initial 
states. Equation (|C5I) in Appendix [Cl provides the contri- 
bution from the quark-gluon scattering subprocess with 
polarized initial-states: 



A'i 



,,19 _ 



8iT 2 a s 
3 



91 

Ql 



1- Wi 



[4 - (z lZ2 - 1) 2 ] 



1 Q 2 

Tq 2 



- r z 2 + (zi + z 2 f] 



x — — 6(s - 



t - 



Z\Z 2 



Aw q L 9 



8ir 2 a s 



Z2 



i + Ql/Q 2 



x T R [z 2 - { Zl + z 2 f 



S 



Z\Z 2 



S(s + t + u - 



A< A 



8TT 2 a s 



i-z 2 \ n 



><2 T « [zl-{zx + z 2 f] 
° -6(s + t + u-Q 2 ) 



Z\Z 2 



Aw q l 



- Aw, 
2 

, 8n 2 a s 



3 

xTji 
S 



Q_ 

Q± 



1- Wl + QW 



ZlZ 2 



■5(s + t + u-Q 2 ) 



(20) 



As shown in Eq. (|C8[) of Appendix [Cl the contribution 
from the gluon-quark scattering subprocess with polar- 
ized initial-states is 



,<i<i 



8TT 2 a s 
3 



[4 ~ i** i) 2 ] 



s 



Z\Z 2 



6(s + i + u-Q 2 ), 



Awf 



8ir 2 a s 



1- Wl + Q^/Q 2 



' 9 3 

x T R [z 2 - {zx + z 2 ) 2 } 



X ■ 



Z\Z 2 



■ 6(s + i + u - Q 2 ) , 



aw AA ~ e g 



87r 2 a s 



i-Wi + qI/q 5 



x-T fl [z\ - ( Zl + z 2 ) 2 } 



S 

z\z 2 



6(s + t + u-Q 2 ) 



— A'i 



Aw<£ 



2 L 
9 87r 2 a s 



9 3 

x T R [z 2 ] 



(^) [i-^i + Qi/Q 5 



S_ 

Z\Z 2 



6(s + t + u-Q 2 ). 



(21) 



We note that other than for Aw A , the contributions 
from the gluon-quark subprocess are effectively the same 
as those from the quark-gluon subprocess, with Z\ and 
z 2 switched. 



C. Limit of Q ± /Q -> 

In this subsection, we examine the analytic behav- 
ior of each parton-level helicity structure function as 
Q±/Q 0. Keeping up to the leading power terms, 
we can simplify the parton-level Mandelstam variables 
and the phase space (S-function as 



Ql 



{l-z x ){l-z 2 ) 

Ql 



Ql 
(i - *i) ■ 



(22) 



The expression for s is an immediate consequence of the 
phase space (5-function, which, in turn, can be expanded 
as [H 

C -Sis + t + u-Q 2 ) 



Z\Z 2 



X\X 2 



6(1 -z 2 ) 6(1 - Zl ) 



(l"*l)+ (1-^2) + 

o 2 

+5(1 - zi) 6(1 - zi) ln^- 



(23) 
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The standard definition of "+" distribution is 



(1-2)- 



= I dz ^M + /(l) ln(l - g) 



(1-z) 



(24) 

Substituting Eqs. {22]) and ([23]) into Eqs. (JUJ) and 
(If 6p . we obtain the analytic behavior of the perturba- 
tively calculated parton-level helicity structure functions 
as Qx/Q — > 0. For the quark-antiquark annihilation pro- 
cess, these are 



Equation ([19]) allows us to conclude that, at this order, 
the analytic behavior of the quark-antiquark annihilation 
subprocess as Qx/Q — > is independent of whether in- 
coming (anti)quarks are spin-averaged or polarized. The 
contributions to the parton-level helicity structure func- 
tions are given in Eq. ([23]) . 

On the other hand, the polarized contributions from 
quark-gluon scattering subprocess are different from 
those for "spin-averaged" initial parton states. From 
Eq. (QD]), we obtain 



w 



AA 



+P qq ( Zl )S(l - z 2 ) 
+2C F 8{1 - z 1 )5{\ - z 2 ) 



H m ] - 2 



+2C F <5(l-zi)<5(l-z 2 ) 



87r 2 a s 
3a;iX2 



P qq {z 2 )5{l - Zl ) 
-P qq (zi)6(l - z 2 ) 



+2C F S(1- Zl )5(l- z 2 ) 



-C i ,[l + z 1 ]<5(l-z 2 )L (25) 



For the quark-gluon subprocess, the small Qx behavior 



87r 2 a s (Q 



q ix\x 2 \Q\ 
87r 2 a 



'!<! 

y AA 



1 2 87r 2 a s 



— e 

2 9 3xiX2 



P OT (*a)«(l-*i), 
P qg (-z 2 )S(l- Zl ), 
P qg (-z 2 )S(l-zi), 



5 *g^(_«) T;i[1 _ 24]j(1 _ 2l)<26) 



Atof 



J AA 



Aw«f 



2 87T 2 Q a / Q 2 

Zq 3xix 2 l,Q2 



AP ra (z 2 )5(l-^) 



2 87r 2 a s 
9 2>x\x 2 



AP qg {-z 2 )8{\- Zl ) 



1 2 87T 2 q; s 

2 9 3xix 2 



AP qg {-z 2 )5{l-z x ) 
8n 2 a s ( Q 



q 3>x\x 2 \Qx 



-T R 6(1- Zl )}, (29) 



where AP qg (z) is the leading polarized gluon-to-quark 
splitting function 



AP qg (z)=T R [z 2 -(l-z) 2 ] 



(30) 



Similarly, based on Eq. ([21]) . the small Qx behavior of 
the polarized gluon-quark contribution is 



Aw 9 T q 



Awf 



•I'l 

AA 



87r 2 a s / Q 



3x\x 2 \Q± 
87r 2 a 



3xiX2 



T^j AP M («i)5(l-«2) 
AP qg (- Zl )5(l-z 2 ) 



A, 



,.m 



2 9 3x1^2 
87r 2 a s ( Q 



AP qg (- Zl )5(l-z 2 ) 



3xix 2 \Q± 



)[T R S(l-z 2 )] 



(31) 



Clearly, the perturbatively calculated helicity struc- 
ture functions at order of a s and beyond are singular 
as Qx/Q — * 0: Wt and Wa have the power divergences, 
Q 2 IQ\ and Q/Qx, respectively, as well as ln(<3/<3x) di- 
vergences, whereas Wl and Waa show ln(<5/<5j_) diver- 
gences da, eh 13 [3 ■ 



The parton-to-parton splitting functions are 



PM 



C F 



1 



(1-*)- 
T R [z 2 + (l-z) 2 



- 2 S(l~z) 



(27) 
(28) 



With zi and z 2 switched, Eq. ([26]) is also true for the 
gluon-quark subprocess, except for u> A 9 which needs an 
extra overall minus sign. Our results for the form of 
the helicity structure functions for unpolarized incoming 
partons as Qx/Q — > in E qs. (1231) and ([2U]) are consistent 
with those derived in Ref. 



III. ASYMPTOTIC CURRENT CONSERVING 
TENSOR 

In this section, we investigate the possible connection 
between the logarithmic divergences of different helicity 
structure functions, and we show that they have a com- 
mon origin. We observe that the four helicity structure 
functions cannot be independent as Qx = where the 
general tensor decomposition in the virtual photon rest 
frame in Eq. (U) is ill-defined. We construct a new asymp- 
totic hadronic tensor that has the right number of inde- 
pendent scalar functions as Qx — ► by requiring that the 
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singular contribution to the hadronic tensor should sat- 
isfy electromagnetic current conservation to all orders in 
a s . We show explicitly that the Q±/Q — > singular con- 
tributions in Wt, Wl, and Waa are related uniquely to 
the singular contribution of the angular-integrated cross 
section. 

The general arguments in Ref. [26[ show that there 
should be only two independent power-divergent scalar 
functions as Q±/Q — > in the Collins-Soper frame. To 
display the explicit dependence of the hadronic tensor on 
Q±/Q, we rewrite the unit vectors of the Collins-Soper 
frame in Eq. as 




X" 



Vz V Q ' 1 



with uniquely fixed. By expanding the full Drell-Yan 
hadronic tensor in Eq. ^ and using Eq. (|32[) in the limit 
Q±/Q —> 0, we obtain the following form for the singular 
terms of the tensor 127, 30] 



1 



V2 



[n\n v + n fl n ± ) e y 

Asym Q Asym j 



+ 



V2 



w. 



Q 



[n\n v + n»n v x ) e" 



Q 



Asym 



Q_ w A s ym 

Q± 



(33) 



At this point, there are two unspecified divergent scalar 
functions: W 2 Asym oc Q 2 /Ql and W^ sym cx Q/Q ± as 



Q±/Q —> 0. In Eq. p3p . the unit vectors n, n, n± specify 
the center-of-mass frame of the hadron collision, defined 
in Appendix [XJ 

The singular tensor as Q±/Q — > in Eq. (f3"3"|) is not 
current conserving since Waging ^ 0- hi order to resum 
the singular terms of the hadronic tensor to all orders in 
a s , we require a tensor that incorporates all the singular 
terms and also conserves the current perturbatively at 
any order of a s . We use the term asymptotic tensor for 
this current-conserving tensor. We define it to be 



Asym 



{-g" v + n»n y + n^n v ) W^ ym 



-9 

Q± 



l 

X 2 



W, 



Asym 



"± 1 o-" 

Q TT/Asym 



Q-i- ( u v i a v i Q-L u i 

-— — n,j_n + wrriy + -prrn rn 



Asym 



-W{ 



Asym 



(34) 



where the components of the virtual photon momentum 
Q + = q ■ n and Q~ = q ■ n are defined in Appendix [A] 
The asymptotic tensor in Eq. (f3"4"|) is equal to the singular 
tensor in Eq. (|33|) plus a minimal non-singular term such 
that^WT =0. 

The angular-integrated cross section is obtained from 
the trace, da/d 4 q oc —g^ v W^. The trace of the asymp- 
totic tensor in Eq. ((34j) should therefore be fixed by the 
asymptotic term W Asynl of the angular-integrated Drell- 
Yan transverse momentum distribution [25j ] . This state- 
ment allows us to fix uniquely the asymptotically diver- 
gent function W^ syra in Eq. J31J). We obtain 



(35) 



The angular-integrated cross section fixes the value of 
W^ sym , but it cannot fix the second scalar function 
^rAsym j n (jg^j) This second function represents the 
singular perturbative behavior of the structure function 
Wa- We defer discussion of Wa until Sec. V and con- 
centrate on transverse momentum resummation for the 
other three helicity structure functions, Wt, Wl, and 
Waa- 

We reexpress the asymptotic tensor in terms of the 
previously defined unit vectors in the Collins-Soper frame 
as 



Asym 



1 

~\ + Q 2 JQ 2 



1 + 



Ql/Q 2 
"l/Q 2 



x»x v 



''XI 

W Asym 



1 



1 



JQ 2 



\X^Z V + Z»X V \ W^ sym . (36) 



Upon comparison with Eq. ([4]), we immediately derive 
the corresponding asymptotic helicity structure func- 
tions, 



W, 



Asym 



l 



l 



JQ 2 



w 



Asym 



w, 



Asym 
AA 



2 1-1 

Ql/Q 2 
i + Q 2 JQ 2 

1 Ql/Q 2 

2 l + Ql/Q 2 



Ql/Q 2 
w Asyn 



|^Asym yyAsyn 



2 

^yAsym 



2 2 

Q2_ M^ Asym 
Q 2 ' 
^ 1 Ql W Asym 
~ 2 Q2" 2 



(37) 



Equation (|37p shows that current conservation relates 
the Q±/Q — > divergent terms of the transverse, longitu- 
dinal, and double-spin flip structure functions intimately 
to the divergent part of the angular-integrated trans- 
verse momentum distribution. The next key question, 
addressed affirmatively in the next section, is whether 
the asymptotic helicity structure functions in Eq. (|3T[> . 
as derived here, are sufficient to remove all the lead- 
ing divergences in the perturbatively calculated structure 
functions order by order in a s . 
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IV. PERTURB ATI VE FINITE TENSOR 

We show in this section that the three asymptotic he- 
licity structure functions presented in the last section in- 
clude all the Q±/Q — > leading divergent terms of the 
corresponding perturbatively calculated helicity struc- 
ture functions, and therefore, that we can define a per- 
turbatively finite tensor from the difference 



W£Y . = WH V 

vr Finite — ''Pert 



wr 

Asym ' 



(38) 



at any order of a s . This finite tensor conserves the cur- 
rent since the asymptotic tensor conserves the current. 

The Q±/Q — > divergent part of the angular- 
integrated cross section is obtained from the trace of the 
hadronic tensor g^ W^ v . Applying this statement at the 
parton level, we use the results of Sec. Ill to derive the 
Q±/Q —> asymptotic terms for the angular- integrated 
and spin-averaged qq — > j*g and qg — > "f*q subprocesses. 
These are 



Asym 
It) - 
w qq 



8TT 2 a s Q 2 



+P qq (z 1 )8(l-z 2 ) 



Asym 

qg 
2 

Asym 

gq 



+2C F S(1 - zi)8(l 



*2 



ln(- 



3x\x 2 Q 



P qg (z 2 )d(l- Zl ) 



2 87T 2 Q S Q 2 

e„- — T P qg (zi)S(l 



" q 3xix 2 Q 



Z2 



(39) 



_L 



Using Eq. (|37p at the parton level, we find that as 
Q±/Q — * 0, the parton-level asymptotic terms in Eq. (f39|) 
remove all divergent contributions of the correspond- 
ing perturbatively calculated helicity structure functions. 
For the quark-antiquark annihilation subprocess, 



1- 



i Ql/Q 2 



2 f + Ql/Q 2 

Ql/Q 2 



Asvm 
W - 
<?<? 



O(Q ± ) 



Asvm 
W - 



Ql/Q 2 
Ql/Q 2 



l 



2 I + Ql/Q 2 



Asym 
II) - 



0(Ql) 
=> 0(Ql). 



(40) 



With z\ and z 2 interchanged, Eq. (|4T|) is also true for the 
gluon-quark subprocess. Other than the non-logarithmic 
finite piece (as Qj_/Q — > 0) in the quark-gluon contribu- 
tions to Wl and VFaa, the asymptotic tensor completely 
removes the leading term of the perturbatively calculated 
helicity structure functions as Q±/Q —* 0. 

The parton-level asymptotic terms for the polarized 
quark-antiquark, quark-gluon, and gluon-quark subpro- 



cesses are 



Aw Asym 



Aw^ ym 



Asvm 

in - 

LU qq 



2 87r 2 a a Q 2 
^3^- 2 Ql AP ^ Z2 

2 8TT 2 a s Q 2 

^3x^x- 2 Ql APq9iZl)5{1 



0*(i-«i). 



- 3») ■ (42) 



Since Ao;^ 



and Au^ sym 



Asvm 

HI - 
1X1 qq 



qq ~ ~qq , Eq. VB 

true also for the polarized quark-antiquark subprocess. 

The finite contributions in the parton-level helicity 
structure functions for polarized quark-gluon or gluon- 
quark subprocesses are not the same as those for the 
corresponding unpolarized subprocesses. We find 



Aw q T 9 



1 



Ql/Q 2 

2 1 + Ql/Q- 



Aw^ ym 



A, 



'■S'-i 



0{Ql) 

Ql/Q 2 
i + Ql/Q 

9 87r 2 a 



A< sy " 
2 

8(1 -zx) 



Ain qg 
AA 



q 3XiX 2 

x[AP qg (-z 2 )-AP qg (z 2 )] + 0(Ql) 
■ ' "I 2 A?u As y m 

w qg 



1 _Ql/Q 2 

2 1 
1 



rjQ 2 



=>■ - e 



S7r 2 a s 



8(1 - Zl ) 



x 



2 q 3xix 2 
[AP qg (-z 2 ) - AP qg (z 2 )] + O(0i).(43) 

With z\ and z 2 interchanged, Eq. ([43]) is also true for 
gluon-quark subprocess. 

The uncanceled finite term in the helicity structure 
functions Wl and Waa is proportional to 



For the quark-gluon subprocess, 



,,'i'J 



,,'i'J 



w 



AA 



1 

1 

2H 
Ql/Q 2 



Ql/Q 2 



_/Q 2 



Asym 

qg 



Asym 

w qg 



i Ql/Q 



w, 



Asym 

qg 



2 I + Ql/Q 



0(Ql 



87r 2 a, 



1 + Ql/Q 2 2 
x[P qg (-z 2 )-P qg (z 2 )} 
i2 irii 



8(1- zi) 



q 3a;ia;2 
-0(Q\) 

> l -e 2 ^8(l 
2 q ixix 2 



x[P qg (-z 2 )-P qg (z 2 )}+0(Ql). 



(41) 



P qg (-z 2 )-P qg (z 2 )=Az 2 T R , (44) 
for unpolarized initial partonic states, and to 



AP qg (-z 2 ) - AP qg (z 2 ) 



-1z 2 T R , 



(45) 



for the polarized initial partonic states. Therefore, for the 
scattering of two polarized hadrons with the same helicity 
(both positive or negative) , the quark-gluon contribution 
to the perturbatively finite term of the helicity structure 
functions Wl and Waa vanishes as Q±/Q — > 0. This 
result is obtained because the perturbative contribution 
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to the longitudinal and double-spin flip helicity structure 
functions is proportional to P qg {— z 2 ) + z "i) m the 

limit of Q±/Q — > 0, the corresponding asymptotic term 
is proportional to P qg (z 2 ) + AP qg (z 2 ), and the difference 
vanishes due to Eqs. (|44)l and (|45|) . We also observe 
that, at this order, the uncanceled term in the quark- 
gluon subprocess is proportional to the helicity flipping 
splitting function, 



P q - g+ {z)=P q+g -{z)=T R {l-zf 



(46) 



The finite term as Q±/Q — > for the quark-antiquark 
subprocess at this order is removed completely by the 
asymptotic term since the helicity flipping splitting func- 
tion for the quark vanishes at this order, P q - q +(z) — 
P q+q -(z)=0. 

Our observations allow us to claim that transverse mo- 
mentum dependent factorization for the full hadronic 
tensor, which is the basis for the Collins-Soper-Sterman 
6-space resummation, breaks at subleading power in the 
Q±/Q expansion, but only in the helicity flipping chan- 
nel. The breaking seems not to supply leading logarith- 
mic terms. 

The asymptotic current-conserving tensor introduced 
in last section is sufficient to remove all leading diver- 
gent terms in the perturbatively calculated hadronic ten- 
sor. The logarithmic terms in the perturbatively calcu- 
lated helicity structure functions, Wt, Wl and Waa, 
are shown here to have the same origin as those in the 
angular-integrated cross section. Therefore, for these he- 
licity structure functions we can obtain a perturbatively 
finite difference as 



VF-Finitc = W Fcrt 



Asym 



(47) 



with i = T,L, AA. 



V. FULL HADRONIC TENSOR INCLUDING 
TRANSVERSE MOMENTUM RESUMMATION 

In this section we present expressions for the trans- 
verse momentum dependence of the structure functions 
incorporating resummation to all orders in a s of the sin- 
gular divergent behavior as Q±/Q — > and including the 
contributions at order a s that are finite in the small Q± 
limit. We begin first with a brief summary of the resum- 
mation formalism developed for the angular-integrated 
cross section. 

As explained above, when Q± <C Q, the Q± distri- 
bution of the helicity structure functions calculated in 
conventional fixed-order perturbation theory receives a 
large logarithmic term, ln(Q/Q±), at every power of a s , 
which is a direct consequence of the emission of soft and 
collinear gluons from the incident partons. Therefore, 
when Q±/Q is sufficiently small, the convergence of the 
conventional perturbative expansion in powers of a s is 
impaired, and the logarithmic terms must be resummed. 



Resummation of the large logarithmic terms can be 
carried out either in Q^-space directly, or in the im- 
pact parameter, 6-space, which is the Fourier conjugate 
of (J_i_-space. It was first shown by Dokshitzer, Diakonov, 
and Troian that in the double leading logarithm approx- 
imation, the dominant contributions in the small Qt re- 
gion can be resummed into a Sudakov form factor [22| . 
By imposing transverse momentum conservation with- 
out assuming strong ordering in the transverse momenta 
of radiated gluons, Parisi and Petronzio introduced a 6- 
space resummation method which allows one to resum 
some subleading logarithmic terms (23|. Using a renor- 
malization group equation technique, Collins and Soper 
improved 6-space resummation to resum all terms as sin- 
gular as ln m (Q 2 /Q\)/Q\, as Q ± -f ||. Using this 
renormalization group improved 6-space resummation, 
Collins, Soper, and Sterman (CSS) derived a formalism 
for the transverse momentum distributions of vector bo- 
son production in hadronic collisions [25[ . This CSS for- 
malism, developed originally for angular-integrated vec- 
tor boson production, casts the cross section in the fol- 
lowing generic form [25| 



da 



1 



(2tt) 2 
Y(Q_l,Q, xi,X2) 



d 2 be i ^ ± ^W(b,Q,x 1 ,x 2 ) 



(48) 



The function W provides the dominant contribution 
when Q± Q, while the Y term supplies contributions 
that are negligible for small Q± but become important 
in practice when Q± ~ Q. The function W in Eq. |48|) 
incorporates all powers of large logarithmic contributions 
from ln(l/6 2 ) to ln(Q 2 ). It has the following form [25[ 



W(b, Q,xi,x 2 ) 



- a S(b,Q) 



W(b 7 c/b, xi,x 2 ) . 



(49) 



where c is a constant of order one 12511, and 



SU,())=I ^4 

c 2 /b 2 M 



In 



^).l(n.J / ,)l+ii(n.J / ,)l 



Functions A(a s ) and B(a s ) may be calculated perturba- 
tively in powers of a s [25]. Function W(b,c/b,XA,XB) 
in Eq. (|4*9"|) depends only on one momentum scale, 1/6, 
and it may be calculated perturbatively as long as 1/6 
is large enough. The large logarithms from ln(c 2 /6 2 ) to 
ln(Q 2 ) in W(b,Q,xi,x 2 ) are completely resummed into 
the exponential factor exp[— S(b, Q)]. The finite Y term 
is defined to be the difference between the cross section 
calculated in conventional fixed-order perturbation the- 
ory and the asymptotic cross section which is equal to 
the perturbative expansion of the resummed part of the 
cross section, the first term on the RHS of Eq. (|4"8")) . 

The function W{b,Q,xi,x 2 ) of the CSS 6-space re- 
summation formalism in Eq. (|48[) is not exactly equal to 
the Fourier transform of the transverse momentum dis- 
tribution, but its Fourier transform reproduces all lead- 
ing divergences of the type ln m (Q 2 /(5 2 L )/(5 2 L in the per- 
turbatively calculated transverse momentum spectrum 
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when Qx/Q —> 0. Combined with the perturbatively 
finite Y term, the Fourier transform of the resummcd 
W(b, Q, x\, Xq) gives a good description of heavy vector 
boson production at collider energies [U, HH • 

The transverse momentum dependence of the angu- 
lar distribution of lcptons from the Drell-Yan mechanism 
is determined by the transverse momentum dependence 
of the helicity structure functions. Only the transverse 
structure function Wt has a leading divergence of the 
type \n m (Q 2 /Q 2 L )/Q 2 L as Q±/Q -> 0. It might be natu- 
ral to consider the resummation of these large logarithms 
into Wt EH, [HI- However, as we demonstrate in 
Eqs. ([51)1 and J2H|), electromagnetic current conservation 
requires that the leading logarithmic divergences of the 
structure functions Wl and Waa share the same origin 
as those in Wt and those in the angular-integrated cross 
section. All are included in one asymptotic function, 
W Asym . Resummation of the large logarithmic terms of 
the Drell-Yan helicity structure functions can therefore 
be accomplished in terms of the resummed contribution 
to the angular-integrated Drell-Yan cross section. Refer- 
ring to Eq. (|37p . we obtain the resummed contribution 
to the helicity structure functions in the Collins-Soper 
frame as 



wz^ - ; 



i Ql/Q 2 w Asym 



' AA 



2 1 



Ql/Q 2 



(53) 



with the asymptotic term in these expressions equal to 
the perturbative expansion of the resummed contribution 
in powers of a s . As in the case of the angular- integrated 
cross section, these expressions are applicable for Q -> 
Qj_. Other effects must be considered when Q± 3> Q [33], 

Substituting the expressions for the helicity structure 
functions Wt and Wl from Eq. ([53")) into Eq. (J2), we 
obtain 



da 



l2ir 3 S 2 Q 2 
+ (2W^ crt 



W 



Rcsum 



w, 



Pert"\ 



(54) 



Using Eq. 
term is 

Y 



, we find that the perturbatively finite Y- 



Utt 3 S 2 Q 2 



(2W| Glt + Wl clt ) - W Asyr 



(55) 



^rResum 



1 



i Ql/Q 2 



W 



Rcsum 



Rcsum 



2 1 + Q 2 JQ 
Ql/Q 2 FT Rcsum 



W 



Rcsum 



AA 



i + Ql/Q 2 

1 Ql/Q 2 

2 1 + Ql/Q 2 



2 ' 

Rcsum 



(51) 



All depend on the same QCD resummed expression 
W-Resum that pertains to the angular-integrated Drell- 
Yan cross section [25]. By comparing Eq. (|4"5)l with 
Eq. (J2J), we obtain 



l2ir 3 S 2 Q- 



W 



Rcsum 



l 



{2-Kf 



d 2 be l Q^W{b,Q,x 1 ,x 2 ). 

(52) 



In analogy to the CSS result for the angular-integrated 
cross section in Eq. ([48]) . the expressions for the full trans- 
verse momentum distribution of the helicity structure 
functions are 



W T = 



W L = 



1 

1 - - 



Ql/Q 2 



2 1 + Ql/Q 2 

Ql/Q 



W 



Rcsum 



Pert / 1 i Ql/Q 2 



2 1 

TxrResum 



i + Ql/Q 2 



2 

Ql/Q 2 



JQ 2 



w Asyr 



W Asyi) 



Waa 



1 + Ql/Q 2 
Ql IQ 2 w Rcsum 



2 1 + Ql/Q 2 



A. Lam- Tung relation 

The Lam- Tung relation states that the longitudinal 
and the double-spin-flip structure functions obey the 
equality W L = 2W A a- Based on Eqs. (37]) and (JSTJ) and 
the definition in Eq. ([53"]) . we find that possible viola- 
tion of the relation can come only from the non-singular 
finite piece of the perturbative contribution. The re- 
summed contribution is known to dominate the angular- 
integrated cross section in the region of small and modest 
Q±, and, by extension, we expect it to dominate the be- 
havior of Wl and Waa in the same region. We conclude 
that violation of the Lam- Tung relation as a function of 
Q± should be relatively small, consistent with the results 
of perturbative calculations at order a 2 [9( , but demon- 
strated here to all orders in a s . 

An alternative way to state the Lam- Tung relation is 
in terms of the angular coefficients A and v, defined in 
Eq. (|A15[) . It is expressed as 1 — A — 2v = 0. We derive 



A = 



Wt-Wl 



pyRcsum 



|^Resum 



W T 
1 - 



W L 



Ty^/" Rcsum Tyj^Resum 



AIQ 2 



i + |Qi/Q 2 

2Waa _ 
Wt + Wl^ 

Ql/Q 2 
i + IQI/Q 2 



2 ^Rcsum 
Tyj/Resum _|_ Tyj^Resum 



(56) 



The analytic expressions in Eq. (|56J) were derived first in 
Ref. 1 based on the perturbative calculation of qq — > 
7*<?. Our result is valid for all orders in a s if we retain 
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FIG. 4: The transverse momentum dependence of the 
angular-integrated Drell-Yan cross section, obtained from the 
contributions of the helicity structure functions, Wt and Wl , 
in Eq. 1|53[) is shown as a solid line and compared with data 
from Fermilab experiment E772 [3t3 ] for Q in the interval 
(8, 9) GeV. The dashed and dot-dashed curves show our cal- 
culations for the contributions from Wt and Wl- The inset 
shows the Wl contribution on an expanded scale. 



only the leading resummed contribution, and it is inde- 
pendent of the type of incident hadrons. 

A recent analysis of Fermilab data shows reasonable 
agreement with the Lam- Tung relation for moderate val- 
ues of Q± [lj|, while earl y d ata with pion beams show 
some violation [H [H, [T3,lll ■ 



B. Phenomenological example 

As an example, we show in Fig. [4] an explicit numerical 
evaluation of the helicity structure functions Wt and Wl 
computed from Eq. (|53p . The double spin- flip structure 
function Waa — Wl/2 since both the resummed contri- 
butions and the finite perturbative contributions at order 
a s satisfy this relationship. We choose the mass interval 
8 GeV < Q < 9 GeV and E hcam = 800 GeV in order to 
compare with data from Fermilab experiment E772 [36[ . 
The parameters we use are identical to those used for 
Fig. 14 in Ref. (H- 

The dashed and dot-dashed lines in Fig. |4] represent 
the Wt and Wl contributions to the cross section, while 
the total contribution is proportional to 2Wt + Wl ■ We 
remark that the transverse momentum distribution after 
resummation is finite as Q± — > for Wt, but it becomes 
vanishingly small in the case of Wl ■ 



C. Discussion of Wa 

As shown in Sec. Ill, the perturbative contribution to 
the single spin-flip structure function Wa is proportional 
to Q/Q±, which is singular as Q±/Q — > 0. Unlike the 
other helicity structure functions, Wa does not show a 
logarithmic divergence in the Collins-Soper frame, a fea- 
ture that seems special for this frame [llj . The absence of 
the divergence could be a consequence of the symmetry 
of the frame with respect to the hadron beam directions, 
which requires W A oc Wie~ 2y - W 2 e 2v in Eq. (|A"9|) . and 
the fact that the leading logarithms arise from the region 
of phase space where z\ — ► 1 and Zi — * 1. 

In this frame, the quark-antiquark contribution to Wa 
is completely antisymmetric in z± and because of the 
opposite sign between the Wi and the W2 terms above. 
The quark-gluon (or gluon-quark) contribution is propor- 
tional to 1 — Z2 (or (1 — z\)). The asymmetry in z\ and z-i 
strongly reduces the numerical size of these contributions 
when Q± ^ 0. 

The combination of the quark-antiquark and 
antiquark-quark subprocesses gives the following 
perturbative contribution to Wa, 

W qq + Wq q (X [q A fa) q B (&) + q A {ti) «b(6)] («i - 4) ■ 

(57) 

This contribution vanishes in the central region for colli- 
sions between hadrons of the same type. The quark-gluon 
contribution also shows a similar asymmetry between z\ 
and z%, 

W qg + W gq (X q A ^l)9B(h) (4 ~ 2zf) 

+.9A(a)to(6)(2z 1 2 -z 2 2 ) . (58) 

Collinear factorization in the perturbative calculation 
ceases to be valid when Q± ~ Aqcd or less. At Q± = 0, 
the helicity structure function Wa itself is ill-defined. We 
might still be able to test the physics of the single spin-flip 
structure function in the small Q± region by introducing 
a new observable, for example, the first moment of the 
structure function, 

W a (Qt,Q)= dQ ± Q ± W A (Q±,Q) (59) 
Jo 

which is perturbatively more stable if Qt is large enough. 

Is it possible that a different kind of resummation 
would handle the non-physical Qj 1 divergence at Q± = 
in Wa? We do not have an answer to this question in 
the collinear QCD factorization approach. However, we 
might gain insight by investigating the angular distri- 
bution from another perspective - starting with trans- 
verse momentum dependent quark-antiquark annihila- 
tion HH3. 

VI. SUMMARY AND DISCUSSION 

Massive virtual photons, the W boson, and the Z bo- 
son have important decay modes into pairs of leptons. 
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The angular distribution of these leptons, measured in 
the rest-frame of the parent states, determines the align- 
ment (polarization) of the massive vector boson and, con- 
sequently, supplies more precise information on the pro- 
duction dynamics than is accessible from the angular- 
integrated rate alone. An understanding of the expected 
angular distribution is also important for estimating cor- 
rections associated with limited angular acceptance in 
typical experiments. The changes expected in the angu- 
lar distribution as a function of the transverse momentum 
Q± of the vector states is a topic of considerable interest, 
both for refined tests of QCD and to reduce systematic 
uncertainties on the determination of the W boson mass 
00. 

In this paper, we calculate the transverse momentum 
Q± dependence of the four helicity structure functions 
for the production of a massive pair of leptons with pair 
invariant mass Q. These structure functions determine 
the angular distribution of the leptons in the pair rest 
frame. We work within the QCD collinear factorization 
approach valid for Q± > Aqcd- Our goal is the pre- 
diction of the full Q± dependence of the four structure 
functions, including the region of small and intermediate 
Q± where the cross section takes on its largest values. 

As also noted by others, when calculated at fixed or- 
der in QCD perturbation theory, the structure functions 
show unphysical inverse-power Qj™ (n = 1 or 2) or loga- 
rithmic \n(Q/Q±) divergences, or both, as Q± — ► 0. For 
the angular- integrated cross section, da/d^q, it is well 
established that similar unphysical divergences can be 
removed after resummation of the ln m (Q 2 / Q\) / Q\ sin- 
gular terms from initial-state gluon emission to all orders 
in a s [I| 11, [II 11. 

We begin our analysis with the observation that the 
four helicity structure functions cannot be independent 
at Q± = 0. The general tensor decomposition in the vir- 
tual photon rest frame in Eq. g]) is ill-defined at Qj_ = 0. 
Then, we employ electromagnetic current conservation 
to construct a new asymptotic hadronic tensor that has 
the right degrees of freedom as Q± — > and embodies 
the minimal divergent behavior present at fixed-order in 
QCD perturbation theory. We find that the leading log- 
arithmic behavior of three of the helicity structure func- 
tions, Wt, Wl, and Waa, has a unique origin. Its ori- 
gin is the same as that of the divergence in the angular- 
integrated cross section. We are able, therefore, to re- 
duce the problem of transverse momentum resummation 
for Wt, Wl, and Waa to the known solution of trans- 
verse momentum resummation for the angular-integrated 
cross section [25[. We prove that the small Q± logarith- 
mic divergences in Wt, Wl, and Waa may be resummed 
to all orders in the strong coupling strength a s , yielding 
well behaved predictions for the Q± dependences that 
satisfy the expected kinematic constraints at small Q±. 
The fourth structure function, Wa, requires a different 
treatment, as discussed in Sec. V.C. 

The main results of our research include the fact that 
electromagnetic current conservation uniquely ties the 



perturbative divergences as Q±/Q — > of the other- 
wise independent helicity structure functions Wt,Wl, 
and Waa to the divergence of the angular-integrated 
cross section. Second, the perturbative divergence in the 
angular-integrated cross section is sufficient to remove 
all leading small Q± divergences of the individual he- 
licity structure functions. Third, transverse momentum 
resummation of the angular-integrated cross section de- 
termines the resummation of the large logarithmic terms 
of the helicity structure functions Wt, Wl, and Waa- 
Finally, the approximate Lam- Tung relation between the 
longitudinal and the double-spin-flip structure functions 
is an all-orders consequence of current conservation for 
the leading perturbatively divergent terms. 

In further work, we intend to examine the Q± de- 
pendence of W and Z boson production, where parity 
violating terms introduce additional helicity structure 
functions. Decay of these intermediate bosons into their 
dilepton channels supplies accurate measurements of the 
masses of the bosons. For W production, more accu- 
rate predictions for the angular distribution of the sin- 
gle observed lepton should complement the missing en- 
ergy technique and lead to an improved determination of 
the mass. The mass of the W boson provides an elec- 
troweak observable that bounds the mass of the Higgs 
boson within the framework of the standard model of 
particle physics [381 ]. 

The use of current conservation to establish connec- 
tions between the divergences of different helicity func- 
tions at Q± — ► in the Drcll-Yan process may have 
immediate application for improving QCD resummation 
and predictions for particle production or other observ- 
ables in semi-inclusive deep- inelastic scattering (SIDIS). 
Unlike the Drell-Yan process, the lepton angles in SIDIS 
cannot be integrated over fully because the measurement 
of the DIS kinematic variables x b and Q 2 requires speci- 
fication of the production angle of the lepton in the final 
state. Like the Drell-Yan cross section, the different he- 
licity structure functions in SIDIS have a ln m (Q 2 /q 2 _) 
perturbative divergence at small values of the particle 
transverse momentum q±, defined in the frame where 
the vector boson and the colliding hadron are aligned 
with each other. All helicity structure functions con- 
tribute to particle production in SIDIS. Only the leading 
singular ln m (Q 2 / q±) / q\ logarithms are resummed in ex- 
isting QCD calculations [2£j, H . Inclusion of the effects 
of resummation for the individual structure functions, 
as described in this paper, should lead to more accurate 
predictions for SIDIS observables, such as particle energy 
flow and rapidity dependence, that could be sensitive to 
the relative size of the different helicity structure func- 
tions. 
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APPENDIX A: DRELL-YAN CROSS SECTION 
AND ANGULAR DISTRIBUTION 

In this appendix, we summarize the basic formalism for 
calculating the cross section for dilepton production in 
the Drell-Yan model and the angular distribution of the 
leptons. The expressions in this appendix also establish 
our notation. 

We consider the scattering of two hadrons of momen- 
tum Pi and Pi , respectively, that produces a virtual pho- 
ton of four-momentum q, A{P\) + B(P 2 ) — ► J*(q) + X, 
that in turn decays into a pair of leptons of momentum 
I and I, as sketched in Fig. [TJ The cross section for this 
Drell-Yan production process can be expressed as 

Jq~dn = 2{2^fsW L ^ wfir (A1) 

The leptonic tensor is 
and the hadronic tensor is defined as 

WV = sJ2(p 1 p 2 \4(o)\x)(x\j v (o)\p 1 p 2 



(A2) 



x 



x(2Tr) 4 S 4 (P 1 +P 2 -q 



= s 



d 4 ze^ z (PiP 2 |4(0) J,(z)|P 1 P 2 ),(A3) 



where 3 ^ is the electromagnetic current. Electromagnetic 
current conservation, g M = 0, and the fact that elec- 
tromagnetic and strong interactions are invariant under 
the parity and time-reversal transformation, allows us 
to express the Lorentz tensor, W^ v , in terms of four in- 
dependent Lorentz scalar functions We choose the 
following four frame-independent scalar functions, 



P^PfWi + P^P^W 2 

i [pf P 2 y + Pf Pf 1 W 3 - gT W 4 . (A4) 



The dimensionless current-conserving tensor and the vec- 
tors are defined as 



fj 



r 2 



fJ>L* 



qfq v 



w P lv j^fS, 



with q^^ v = 0. Our choice of the four frame- 
independent scalar functions is slightly different from 
that in Ref . [3] ■ We find that this choice is convenient for 
connecting to the parton-level perturbative calculation 
discussed below. 

By contracting the leptonic tensor L^ v and hadronic 
tensor W M „ in Eq. (|A1[) , we can express the Drell-Yan 
cross section in terms of the four scalar functions Wj and 
the measured hadron and lepton momenta. 

The physical meaning of the scalar functions can be ap- 
preciated if we express them in terms of the four indepen- 
dent "helicity" structure functions, Wi with i — T,L, A, 
and A A, corresponding to the transverse spin, longitu- 
dinal spin, single spin flip, and double spin flip contri- 
butions to the Drell-Yan cross section [J. The helicity 
structure functions are defined in the dilepton center-of- 
mass frame (the virtual photon's rest frame). 

The full hadronic tensor in Eq. (|A4[) can be also writ- 
ten in terms of the helicity structure functions and unit 
vectors in the virtual photon rest frame as in Eq. ^ [i[ . 
In this frame, the lepton momenta are 



I" 



Q 



(1, sin 6 cos cf>, sin 6 sin 0, cos( 



= — (1, — sin6*cos0, — sin^sin^, — cos^) . (A6) 

Substituting the hadronic tensor in Eq. Q and the lep- 
tonic tensor in Eq. (|A6|) into Eq. (|Aip . one gets the dif- 
ferential cross section of Eq. ([I]) . 

The frame- independent structure functions and the he- 
licity structure functions are uniquely related to each 
other once we make a choice of the coordinate system, or 
the unit vectors, in the virtual photon rest frame. The 
unit vectors for the Collins-Soper frame are chosen as (2(| 



Z 1 ' 

X 1 * 
yf 



VQ 2 + Q'i 
Q_\ 

e^T u Z a X . 



q P2 Pf - q Pl Pf 



q P2 Pf + q Pl P£ 



(A7) 



The dimensionless current-conserving hadron momenta, 
Pf and Pf , are defined in Eq. |[A"5]) . and q Pi = P r q/^/S 
with i = 1,2. The hadron and the virtual photon mo- 
menta can be expressed in the center-of-mass frame of 
the collision as 



pf _ 
r 2 ~ 



n' 



Qj 



n 



± ■ 



(A8) 



with total center-of-mass collision energy 
^{Q 2 + Ql)/2 e y, and Q~ = ^/(Q 2 + Q\)/2e~y. In 



Eq. 



= n» = 5^-, and n>l = S^ 1 - are unit 



vectors that specify the light-cone coordinates of the col- 
lision center-of-mass frame, with n 2 = n 2 = 0, n\ = — 1, 
(A5) n ■ n — 1, and n± • n = n± ■ n = 0. In the Collins-Soper 
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frame, the helicity structure functions can be expressed 
in terms of the frame-independent structure functions in 
Eq. ((All) as 



W T = Wa + 1 



2 Q 2 



- (Wi c- 2y + W 2 e +2y ) + -W 3 



W L 



- (Wi c- 2y + W 2 c +2y ) - Jw 3 + W 4 . 



Waa = 



-(me 



W 2 e 



+2y\ 



(A9) 



with the coefficients of the angular dependence given by 

W T - W L 



A 



W T + W L ' 

W A 
Wt + Wl 1 

2IUaa 
Wt + W l ' 



(A15) 



APPENDIX B: PERTURB ATIVE 
CONTRIBUTIONS FROM UNPOLARIZED 
PARTONIC STATES 



From the QCD collinear factorization formalism for the 
hadronic tensor in Eq. we obtain similar factorizcd 
relations for structure functions, 

ab 



with i = T,L, AA, A; and 

ib 



6 



(A10) 

0a(&)&(6)£*Kl,6,«) (Ail) 



with i — 1,2,3,4. Using Eq. (|A9|) . we derive the corre- 
sponding relation between the short-distance parton-level 
structure functions: 



Wt 

wl 

WAA 
WA 



0->4 



IQI 
2 Q 2 



l(u; 1 e~ 2y +u J2 e +2y ) + ^ 3 



- (uji e 2y + oj 2 c +2y ) - -co 3 + uj 4 . 



2 Q 2 



Q± 

Q 



i 



l 
4 

• uj\ e 



(wi e~ 2y 



+ oj 2 e +2y ) + -uj 3 



-2y 



■ uj 2 e 



-2y 



(A12) 



Integration over the solid angle of the decay leptons 
gives the angular-integrated Drell-Yan cross section, 



da 



l2ir 3 S 2 Q< 



12w 3 S 2 Q 2 



(2W t + W l ) 



(A13) 



One can write the normalized Drell-Yan angular distri- 
bution as 



dN_ 

~dn 



da 



da 



\<Pq ) d'qdO 
3 / 1 



4tt V A + 3 



1 + A cos^ 



(A14) 



+(i sin(20) cos <p + — sin 9 cos(2^>) 



In this appendix we summarize the perturbative con- 
tributions to the parton-level helicity structure functions 
for unpolarized initial-state partons. 

Using the definition in Eq. , we derive the contribu- 
tion to the parton-level hadronic tensor from the quark- 
antiquark annihilation diagrams in Fig. [2] with unpolar- 
ized initial parton states. 



M „ 4 2 8n 2 a s 



4£ 2 Q 2 S P^P{ - A&Q 2 S P^Pz 



-{{Q 2 ~t) 2 + {Q 2 -u) 2 )rr 

x S5(s + t + u-Q 2 )i 



(Bl) 



where 4/9 = (l/3) 2 ^A Tr [* 1 \ is the color factor with 
SU(3) generator t A , and 87r 2 a s = (2n)g 2 . The factor 
(2ir) comes from the phase space expression 



S(2tt) 4 6 4 ( Pi 



p 2 -q-pi) 



d 3 P4 



(2n) 3 2E 4 



= 2ir S S(s + i + u - Q 2 ) . 
The parton-level Mandelstam variables are 

s = (pi +p 2 ) 2 = Ci6 s , 

t = (pi-q) 2 =Q 2 -2£iPi-q, 
U = (pa - qf = Q 2 - 2&P 2 ■ q . 



(B2) 



(B3) 



Using Eqs. (jBip and (| Al 2(1 . we obtain the parton-level 
frame-independent structure functions 



,n - 



'3 



4 2 8ir 2 a s 
9 E * 



A -e 2 
9 9 



tu 
8ir 2 a s 
til, 



-4£ 2 Q 2 S] S5(s + t + u-Q 2 ), 
-4£ 2 Q 2 S] SS(s + i + u-Q 2 ). 



4 2 8ir 2 a s 
9 6 " 



tu 



2 e -2y 



e 2 e 2y ] (Q 2 + Qi)S 



xSS(s + i+u-Q 2 ) 



(B4) 



and the corresponding parton-level helicity structure 
functions in the Collins-Soper frame, 



W 



qq 



8-K 2 a s 



til 



[eie- 2y +e 2 e 2y ] S(Q 2 + \q\) 



10 



w 



,fli — 



.</</ 



AA 



W ■ 



xSSis + t + u-Q 2 ), 

y tu 

xS6(s + t + u-Q 2 ), 

x S5{s + i + u-Q 2 ) 

o W L > 



4 e?^ [-^e-^+Cle 2 "] (SQ 2 ) 



ft) 



x S5(s + £ + u-Q 2 ) 



(B5) 



From the quark-gluon scattering diagrams in Fig. [3] 
with unpolarized initial parton states, we derive the 
quark-gluon contribution to the parton-level hadronic 
tensor 

1 , 8n 2 a s 



6 « S(-t) 



8£ 2 Q 2 S Pf Pf - 4£f Q 2 S P^ 

-46 6 2 5 [Pf PT + PfPfl 
-((0 2 -t) 2 + (g 2 -j) 2 )^ 

x5i(s+f + il-Q 2 ), (B6) 

where 1/6 = (1/3) (1/8) £a Tr[i A f/ 4 ] is the color factor. 
We obtain the parton-level frame-independent structure 
functions 

- i^^l [-8^ 2 Q 2 5] S6(s + t + u-Q 2 ), 



6 « $(-*) 



W, 



1 2 8^ 2 a s ^ q2 ^ s s{§ + i + „ _ q2) 



6 9 s{-t) 

- l^ 8 ^ [_86 6 Q 2 S] S6(s + t + u-Q 2 



6 « £(-£) 

1 2 87I " 2q; s 

wl = -e 



6 9 S(-t) 



£ 2 e"^ 5(Q 2 + Ql) 



-(o 2 -a6s) 2 



55(s + t + u-g 2 ); (B7) 



and the corresponding contribution to the parton-level 
helicity structure functions in the Collins-Soper frame, 



1 2 8n 2 a s 
— ez 



6 9 s(-t) 



e 2 e- 2 ^ Sg 2 + (Q 2 - 66S) 2 

-^fe 2 e 2y + 2a6] 5 Q 2 
xSS(s + i+u-Q 2 ), 



1 , 87r 2 a s 



6 9 s(-t) 



2 p- 2 !/ 4- £ 2 p 2 ^ 



^e^ + 2^2] 



w 



AA 



1 2 87r 2 a s 
- e„ — 

6 9 $(-f) 



x (SQi) 55(s + < + w-g 2 ) 

[2e 2 e- 2 ^+^ 2 e 2 « + 2e 1 e 2 ] 



SQl S5(s + t + u-Q 2 ) 



1 



'/.'/ 
2~ L 



2 W 

ie 2 
6 9 



2 e -2 y 



(sg 2 ) S5(s + i- 



Qi 
i 

g 2 ) • (B8) 



Similarly, we derive the contributions to the parton- 
level hadronic tensor from the gluon-quark scattering di- 
agrams. They are the same as those from the quark-gluon 
scattering diagrams with the momenta p\ and p 2 (or 
equivalently with t and u, and £i and £2) interchanged. 



APPENDIX C: PERTURB ATIVE 
CONTRIBUTIONS FROM POLARIZED 
PARTONIC STATES 

In this appendix we summarize the perturbative con- 
tributions to the parton-level helicity structure functions 
for polarized initial-state partons, defined as the states 
with incoming parton polarization projected onto the dif- 
ference of the parton helicity states. 

Based on the same quark-antiquark annihilation dia- 
grams in Fig. [21 we find at this order that the contribution 
to the parton-level hadronic tensor from the scattering of 
a polarized incoming quark and antiquark is the same as 
that from the scattering of an unpolarized quark and an- 
tiquark, 



LLV 



(CI) 



On the other hand, the quark-gluon scattering dia- 
grams in Fig. [3] with polarized quark and gluon initial 
states give a contribution to the parton-level hadronic 
tensor that differs from that for scattering of an unpolar- 
ized quark and gluon, 



(C2) 



We derive 



1 2 87r 2 a s 
6 69 



s(-t) 



^e 2 Q 2 sp^p^ 

-^i&Q 2 S{P?P 2 » + 
-{{Q 2 -i) 2 -(Q 2 -s) 2 W 
: S6(s + t + u-Q 2 ). (C3) 



The contributions to the parton-level frame-independent 
structure functions are 



Awf = 0, 



Aw, 



Aw? 



1 2 87r 2 a 

1 n 8-K 2 a s 

— e„ — 

6 9 s(-t) 



' r 4£ 2 2 g 2 S| S6(s + t + u-Q*), 



[8£i 6 g 2 S] S5(s + t + u-Q 2 ), 
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1 9 8ir A a 
6 69 



2y 



S(-t) 
(Q 2 - 66S) 2 



S(Q 2 + Qi) 



55(s + t + u-Q 2 ). (C4) 



The corresponding contributions to the parton-level he- 
licity structure functions in the Collins-Soper frame are 



A'i 



,.93 _ 



1 2 8n 2 a s 



e 2 e- 2 ^Q 2 -(Q 2 -^) 2 



+l^kmc- 2 y + e 2 c 2 y + 2^ 2 ]sQ 2 

x S5(s + t + u-Q 2 ), 



A; 



n19 _ 



1 p2 87r 2 a, ^ ^ 



266 



A, 



QQ 
,J AA = 



s e 
i( 



6 9 S(-t) L s2 

x (S*g 2 L )S'J(s + t + u-Q 2 ) 



6 9 s(-l) 



2 L 

6 £ i s(-t) [ ^ 2 1 

x (S*g 2 )S'(5(s + t + u-Q 2 ). 



(C5) 



Similarly, we derive the contribution from the polarized 
gluon and quark scattering process, 



99 6 9 s(-u) 



+ 4£ 2 Q 2 SP?PZ 

+46 6 <2 2 5 [PfP^ + P???} 
-{{Q 2 - uf - (Q 2 - s) 2 W» 
xS5(s + t + u-Q 2 ). (C6) 



The contributions to the parton-level frame-independent 
structure functions are 



Atuf 


= ie 2 
6 9 


87r 2 a s 


s(-u) 




= 0, 




Ac 99 


= ie 2 
6 9 


8ir 2 a s 


Auf 


= ie 2 
6 9 


87r 2 a s 


s(-u) 



[46 2 Q 2 S] S6(s + t + u-Q 2 ) 



[86 6 Q 2 S] S5(s + t + u-Q 2 ), 



6 2 e 2y 5(g 2 + Qi) 

5(5(1 + t + u-Q 2 ). (C7) 



The corresponding contributions to the parton-level he- 
licity structure functions in the Collins-Soper frame are 



Aw 99 



1 , 8n 2 a s 



6 9 



6 2 e 2 ^Q 2 -(Q 2 -66^) 2 



^[ei e " 2y + 2 ^ 2 e 2 ^ + 266] SQ 2 



2 |*l " 1 -, »2 

2^ 



2 Q 

x S^s + i + u-Q 2 ), 



A'i 



,,99 



6 * s(— u) 

x (s'Qijs'JOs + i + fi-Q 2 ). 



A„,99 



1 , 87r 2 a, 



— e 



6 9 s(-fi) 



1 



[-C 2 



2„-2y 



266 



A; 



QQ 



[ - SQ 2 ± ) S6(s + t + u 



Awf 

1 , 87T 2 a, r ^o _9„n / Q 



6 -1 s(-u) ^ G ^ { Q 

x (SQ 2 ) S5(s + t + u-Q 2 ). 



(C8) 
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